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Introduction 


At  present,  in  solving  problems  of  turbulence  arising  in  boundary  layers,  the  study  of  their 
receptivity  to  external  disturbances  acquires  great  significance.  In  this  case,  disturbance  genera¬ 
tion  in  boundary  layers  is  attributed  to  stationary  and  nonstationary  factors  (surface  vibration, 
acoustic  waves,  external  turbulence,  etc.).  One  of  the  first  authors  to  treat  this  problem  and  to 
introduce  the  term  "receptivity"  was  Morkovin  [1].  There  are  a  number  of  experimental  and 
theoretical  studies  of  subsonic  boundary-layer  receptivity.  A  detailed  survey  can  be  found  in  [2, 
3].  Much  less  is  known,  however,  about  the  receptivity  of  a  supersonic  boundary  layer,  with  all 
the  publications  in  this  field  being  devoted  to  studies  of  interaction  of  acoustic  waves  with  the 
boundary  layer  [4-6].  At  the  same  time,  a  supersonic  flow  may  contain  not  only  acoustic  distur¬ 
bances,  but  also  vortical  and  entropy  or  heat  waves.  The  interaction  of  those  with  a  supersonic 
boundary  layer  has  not  yet  been  adequately  investigated.  The  problem  of  receptivity  of  vortical 
disturbances  at  subsonic  velocities  has  been  analyzed  in  some  publications,  but  there  is  no 
agreement  among  the  authors  about  the  mechanism  of  disturbance  generation  in  the  boundary 
layer.  In  this  work,  excitation  of  stationary  and  nonstationary  waves  by  external  hydrodynamic 
disturbances  in  subsonic  and  supersonic  boundary  layers  was  calculated. 

Formulation  of  the  Problem  and  Governing  Equations 

A  linear  statement  of  the  problem  is  considered.  A  compressible  gas  flow  in  the  boundary 
layer  on  a  flat  plate  is  taken  as  the  initial  undisturbed  flow. 

Following  Petrov  [7],  the  disturbances  in  the  boundary  layer  are  considered  in  the 
orthogonal  system  of  coordinates  (i;,\| /,z)  fitted  to  the  stream  surfaces  of  the  main  flow.  Here 

y/  is  the  stream  function,  (^,\| i,z)  for  the  plate,  Re  =  ^Ju^x/  is  the  Reynolds  number. 

Moo  >  voo  are  ^e  velocity  and  kinematic  viscosity  of  the  incident  flow,  and  x,  y,  and  z  are  the 
longitudinal,  normal-to-wall,  and  transversal  coordinates  of  the  Cartesian  coordinate  system 
with  the  origin  on  the  plate  edge.  The  gas  is  perfect  with  a  constant  Prandtl  number  Pr.  Making 
use  of  the  estimates  with  respect  to  the  integer  powers  of  the  Reynolds  number  Re  as  in  [8]  and 

rejecting  terms  of  order  Re  with  respect  to  the  dominant  terms  of  the  linearized  Navier  - 
Stokes  equations  [7]  for  disturbances  of  the  form  a(i;,\|/)exp(/ai;  +  /p:r-/(Dt),  we  obtain 


02V  =-(02  In p)v  - [ia - (0j lnM)+0j ]w  -i$w-ucp I p- gmud\p  +  ud-\^T /t), 


d2\_p  +  2|i(/ai7  +  /piv- 2e0  / 3)]  =  -p(h\u  +  dt )v  +  iax\2  +  ;P^23  > 


02  x12  =  (/a  +  0j )p  +  p (02w)v  +  p(0pi  +  df)u  +  u(d]u) p  - tax) i  - i'P  x13 , 


02m  =  -(/a  +  0j)v  -(02w)p/p+  x12  / p 

02x23  =  i$p  +  pdtw-iaxn  -  /'Pf33 ,  (1) 
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d2w  =  -z'Pv  +  x23  /  p  , 

d2q  =  i(op  +  \p(d2H)~  ia  |i(02z./)]v  +  (dj//)(pz7  +  z/p)+(a2  +p2)|i h  /Pr+pdtH  -u(iaxu  +  i  P^b) 
d2h  =  -  Pr(52z/)z7  -  (52/?)fi  /  p  +  Pr(g  -  z/x12 )/ p , 

where  dx  =  ;  52  =  p u  i  dt  =  wc  +  M^1 5  wc  =  ,aM  ~~  ,ro  ;  /q  =  -d\  ln(pzz) ; 

xn  =  2p(zaz7-e0  /3);  X33  =  p(zaw+z'Pz7);  X33  =2p(/pw-e0/3);  Tq  =  ~(32  lnp)v-Kcp/p  ! 
u  is  the  velocity,  T  is  the  temperature,  p  is  the  density,  H  =  h  +  u  /  2  is  the  total  enthalpy,  p 
is  the  viscosity,  p  is  the  pressure,  v  ,  w  are  the  complex  amplitudes  of  the  transversal  and  nor- 
mal-to-stream-surface  components  of  velocity  disturbance,  p/p  =  gmp~T /T ;  T  =  gm\h  ; 
H  -h  +  uu  ;  gm  =  1  /  p  ;  gml  =  l/cp  ;  ( cp  is  the  specific  heat  capacity  of  the  gas  at  constant 
pressure).  The  form  of  the  equations  is  not  changed  by  normalization  to  the  following  scales: 
Voq  / Moo  f°r  length,  Vgo  / um  for  time,  p^  for  viscosity  and  the  stream  function,  for  veloc- 
ity  and  its  disturbances,  Tm  for  temperature,  p^  for  density,  u(D  for  enthalpy,  p^M^,  for 

pressure  and  disturbances  of  viscous  stresses,  p for  the  quantity  q  ,  and  /  Tm  for  spe¬ 
cific  heat  capacity  (the  subscript  00  refers  to  the  values  in  the  incident  flow).  In  this  case,  we 
have  gffl  =yM2  and  gm\  =(y  — 1)M" ,  where  y  =  cplcv  is  the  ratio  of  specific  heats,  and  M  is 
the  Mach  number. 

We  substitute  the  independent  variables:  Re  =  ^,  dr\  =  elf  / u ,  /  =  \|//Re.  Then,  we 
obtain  d{d  =  (l/Re)(S7  +  fxd'),  d2d  =  pa'/Re,  where  0  =  ^5/5R e  ;  the  prime  stands  for  the 

derivative  with  respect  to  13,  and  /j  =-/ /(2ReM).  In  this  case,  the  equations  based  on  the 
estimates  of  the  critical  layer  have  the  following  form: 

v'  =  -gmuTSp  +  p  Tv  -  T(f0u’  +  5)7  -  uw  -  icTr  -  ( f2pT '  -  ud)f  -  fxTu'  +  f2f' , 
p'  =  ~{’c  +  rhl,)v  +  ix\l  +  4^23  “  2P;.«w  > 

*12  =(ix+T8)p  +  pu'v  +  ( ic  +  f\u'  +  ud)u  +  f2u'r  -  iT  +  f2u ' , 
u'  =  —ixv -u'ptT  +  x12/p;., 

^23  =  ’zP  +  Oc  -  Pa  +  u8)w-iz\irUw  +  f2w'  , 
w'  =  -4v  +  x23/pr,  (2) 

q'  =  icoRTp  +  pHv  +  [icu  +  fH'  +  f2u'  +  ird^t  +  (ic  - p a  / Pr+  ud)h  +  f2H7-  - uit  +  f2uu'  +  f2h ' 

h  '  =  -Pvu'u  -h'pfT  +  Pr(</ -mxi2)/ ps , 

Here  uw  =  ixu  +  izw ;  it  =  ixpruw  +  pai7  ;  pa  =  (zv2  +  i:~ )i,. ;  p  =  n- 2p(/az7  +  z'P>v - / 3) ; 

r  =  p  /  p  =  gmp  -  p T  ;  ic  =  Re  uc  =  i  Re(z/a  -  a>);  ix  =/aRe7 ;  iz  =z'PRer  ; 

rh  =R e/ij  =f0u'  +  fxpT' ;  /0  =  -/]  lu  ;  f2  =  f\u  ;  p;.  =pp/Re;  and  pr  =  d\n\xl dT  . 

With  allowance  for  the  substitutions,  systems  (1)  and  (2)  have  the  form 

Z'  =  [A  +  Dd)Z ,  (3) 


57 


where  Z  =  (p,  v,u,w,h,  x12,  123,(7)  in  system  (2)  and  A,  D  are  the  square  matrices  of  the 
specified  functions  Re  ,  q  . 

The  parabolized  system  (3)  is  solved  under  the  following  boundary  conditions.  The 
disturbances  of  velocities  and  temperature  (or  heat  flux)  on  the  surface  are  zero: 

v(o)  =  m(o)=  w(o)=  7’(o)  =  0  (r(0)  =  0).  (4) 

Outside  the  boundary  layer,  the  disturbances  are  obtained  from  their  corresponding  free- 
stream  values  (without  the  model). 

Numerical  Scheme  and  Results 


On  using  the  approximation  j^e  ~  (a  -  «o)/ARe  ( A  Re  =  Re-Reo  is  the  step  of  the 

marching  scheme  and  the  subscript  0  corresponds  to  the  values  calculated  at  the  previous 
step),  system  (3)  is  transformed  into  the  system  of  ordinary  differential  equations 

Z'  =  AZ  +  B(Z-Z0).  (5) 

Its  general  solution  is  a  superposition  of  four  linearly  independent  solutions  of  a 

4 

homogeneous  system  and  an  arbitrary  quotient:  Z=  ^Cm(.r)Zm  +  Z5 .  Herewith,  the 

m= 1 

coefficients  Cm  are  found  from  the  conditions  on  the  surface  (4).  The  selection  of  four 
fundamental  solutions  should  be  made  from  those  decaying  as  q  — >  00  .  The  quotient  should  be 
a  solution  that  describes  external  (undisturbed  by  the  boundary  layer)  waves.  Four  fundamental 
vectors  were  choose  as  following: 


Zl 


=  (0,-1, 0,0, 0,1 

z4  = 


-;a,-/p,l),  Z2  =(2X,  1,0, 0,0, 0,0, 0),  Z3  =  (0,-l,0,0,0,-i'a,-i'P,-l/gml), 
(°-  \ gml ,  iugm  1 ,  /Pgmi  ,0,2iaXgml  ,2/pXgml  ,2iaXgml ) , 


Explanation  of  such  a  choice  you  can  find  in  [14]. 

In  the  present  paper,  the  case  is  considered  where  the  external  disturbances  are  entrained 
by  the  flow.  In  addition,  the  fifth  vector  Z5  corresponding  to  external  disturbances  should  be 
specified  at  the  boundary-layer  edge.  The  work  was  mainly  devoted  to  the  interaction  between 
a  boundary  layer  and  external  non-vortical  hydrodynamic  disturbances,  in  this  case 
Z5  =(0,— Xgml,iagml,ifigml,0,- 2iaXgml ,-2ifiXgml,-2iaXgml) ,  it  was  normalized  to  the 


fourth  component. 

Figure  1  shows  the  amplitude  distributions  of  longitudinal  velocity  disturbances  over  the 
boundary  layer  for  various  values  of  the  frequency  parameter  at  p  =  10  .  Figure  1  also  shows 


a  dependence  proportional  to  q u' ,  obtained  analytically  in  [9]  for  M  =  0,  which  is  a  self-similar 
eigenfunction  of  the  boundary  layer  equations  for  small  values  of  p  [10].  Similar  results  for 


M  =  0  were  obtained  in  [11,  12]  by  the  method  of  optimum  disturbances  specified  at  1=  X()  in 
the  absence  of  disturbances  at  the  boundary-layer  edge.  Weak  dependence  of  the  profile  of  |z7| 
on  changes  in  p  and  disturbances  for  x  =  X0  is  discussed  in  [12].  The  independence  of  the  am¬ 
plitude  profile  u  on  the  parameters  of  disturbances  at  the  external  edge  of  the  boundary  layer 
seems  to  be  caused  by  factors  mentioned  in  [12].  The  main  of  those  factors  is  the  fact  that  exter¬ 
nal  disturbances  generate  eigen  disturbances  with  an  amplitude  u  ~  yu  . 
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Figure  2a  shows  the  maximum  values  of  the 
longitudinal  velocity  amplitude  inside  the  bound¬ 
ary  layer  versus  the  Reynolds  number  for 

co  =  10~6  ,  M  =  0,  and  various  wavenumbers  P  .  It 
can  be  seen  from  Fig.  1  that  the  results  coincide 

for  co  =  10~6  and  co  =  10~7  .  This  implies  that  the 
disturbances  may  be  considered  stationary  for 

co  <  10~6 .  It  follows  from  the  results  shown  in  Fig. 
2a  that  the  maximum  values  of  the  amplitudes  cor¬ 
respond  to  wavenumbers  (P  Re)!  =0.7  for  fixed 
co  and  (PRe)2  =0.55  for  a  fixed  Reynolds  num¬ 
ber  (wavenumber  based  on  the  thickness 
tc/Re  =  Re  ).  These  data  are  in  good  agreement 
with  the  results  obtained  by  the  method  of  opti¬ 
mum  disturbances  [11]  excited  in  the  region  of  the 
flat-plate  leading  edge.  In  addition,  the  value 
(PRe)2  =0.55  is  in  good  agreement  with  the  ex¬ 
perimental  results  on  evolution  of  streamwise 
structures  [13].  It  was  found  in  [13]  that  the  z  pe¬ 
riod  is  approximately  equal  to  the  doubled  thick¬ 
ness  of  the  boundary  layer,  i.e.,  7.  »  2  ■  5x  /  Re  , 
therefore,  /?  Re  =  In  / 10  *  0.6  (it  is  taken  into  ac¬ 


\u 


1  2  3  4  5  6 

r\ 


Fig.  1.  Distributions  of  the  amplitudes  of 
longitudinal  velocity  disturbances  over  the 
—3  —5 

boundary  layer  for  (3  =  10  andco  =  10 

and  M  =  0  (1),  co  =  l(T6  and  M  =  0(2), 

co  =  10  7  and  M=0  (3),  and  co  =  10  6  and 
M=2  (4);  curve  (5)  is  the  dependence 
proportional  to  r\n'  [9]. 


count  here  that  x  =  Re  ). 

The  calculation  results  for  M  =  2  and  co  =  10~6  are  plotted  in  Fig.  2b.  Qualitatively,  they 
are  close  to  the  results  obtained  for  M  =  0.  However,  the  amplitude  of  disturbances  due  to 
external  hydrodynamic  waves  decreases  if  the  boundary  layer  is  supersonic. 


Fig.  2.  Maximum  values  of  the  longitudinal  velocity  amplitude  inside  the  boundary  layer  versus  the 
Reynolds  number  for  M  =  0  (a)  andM  =  2(b),  to  =  10  6  and  (3  =  0.4-10  7  (1),  (3  =  0.6  -10  7  (2), 
(3  =  0.8  ■  10  7  (3),  (3  =  10“3  (4),  and  |3  =  1.2-10~3  (5). 
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max  |w| 

,4 


Figure  3  shows  the  dependences 
max|7|(Re)  for  (3  =  10  ^ .  It  follows 


0  | - , - 1 - . - 1 1 - ' - ' - ' - 1 

100  300  500  700  900  Re  1100 


from  Fig.  3  that  the  amplitude  signifi¬ 
cantly  depends  on  the  frequency  parame¬ 
ter  co :  as  ©  increases,  the  amplitude  of 
disturbances  decreases. 

From  the  analysis  of  the  phase  ve¬ 
locity  c  versus  Re  (obtained  from  the 
phase  growth  along  the  x  axis  for  the  dis¬ 
turbance  u  in  the  vicinity  of  its  maxi¬ 
mum,  i.e.,  for  r|  =  2.3 )  for  a  =  0.3  - 10  4  , 


Fig.  3.  Maximum  values  of  the  longitudinal  velocity 
amplitude  inside  the  boundary  layer  versus  the  Reynolds 

-3  ,  a-4 


number  for  M  = 


(3  =  10 


and  co  =  10 
6 


(1), 


10  (2),  and  co  =  10  (3) 


P  =  10  3  and  M  =  0  and  2,  it  follows  that 
the  phase  velocity  varies  from  0.7  to  0.9 
over  the  range  of  Re  =  200  4-  400  .  In  ex¬ 
perimental  studies  of  propagation  of  in¬ 
cipient  turbulence  spots  [3],  it  was  found 
that  their  leading  edge,  which  is  in  the 
range  of  high  Re  values,  propagates  with  a  velocity  c  «  0.9  ,  and  the  trailing  edge  propagates  with 
a  velocity  c  ~  0.9  ,  which  is  in  qualitative  agreement  with  the  results  of  our  calculations.  It  should 
be  taken  into  consideration  that  the  phase  velocity  was  calculated,  whereas  the  object  of  calcula¬ 
tions  in  experimental  studies  was  most  probably  the  group  velocity.  The  present  calculations  show, 


however,  that  «  1  -  and  the  group  velocity  is  close  to  the  phase  velocity. 


As  it  has  been  mentioned  above,  the  results  presented  in  this  paper  were  obtained  for 
nonvortical  external  disturbances.  Additional  calculations  allowing  for  vortical  disturbances  in 
the  external  flow  entrained  by  the  main  stream  show  that  the  character  of  disturbances  inside 
the  layer  remains  unchanged  in  this  case.  Only  the  relation  between  the  amplitudes  of  external 
and  internal  disturbances  undergoes  some  changes.  The  results  obtained  can  be  explained  as 
follows.  For  high  Reynolds  numbers,  the  evolution  of  disturbances  does  not  significantly 
depend  on  the  presence  of  waves  at  the  boundary-layer  edge,  and  for  low  Reynolds  numbers, 
the  amplitude  of  disturbances  inside  the  layer  is  determined  by  conditions  at  x  =  Xq  .  A  change 


in  external  conditions  at  the  boundary-layer  edge  is  equivalent  to  a  change  in  the  initial  data  at 
x  =  Xq  ,  which  results  in  variation  of  the  absolute  value  of  the  amplitude  for  x  >  Xq  ,  whereas 


the  character  of  the  dependence  of  |w|  on  rj  is  determined  by  the  first  eigen-solution  [12]. 

It  was  also  researched  in  this  paper  an  interaction  between  a  boundary  layer  and 
temperature  disturbances.  In  this  case  we  took  Zg  =  (0,0,0,0,l,0,0,->./Pr) . 

Some  results  you  can  see  at  the  following  figures: 

Figure  4  shows  the  distribution  of  velocities  and  temperature  over  boundary  layer  (1  - 
500|v| ,  2  -  5|m|,  3  -  500|w|,  4  -  |/7| ;  ^0.5-KT3,  ®  =  1(T6,  M  =  2,  Re  =  103),  it  can  be 

seen  that  amplitude  of  |m|  is  more  than  10  times  lower  than  \h  I .  The  shape  of  |m| -profile 
remains  the  same  as  it  was  in  fig.  1 . 
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Fig.  4.  Fig.  5. 


Figure  5  shows  maximum  values  of  the  longitudinal  velocity  amplitude  inside  the 
boundary  layer  versus  the  Reynolds  number  for  temperature  disturbances  case.  M  =  2, 

®  =  10“6  and  /?  =  0.4-l(T3  (1),  /?  =  0.6 -10“3  (2),  y?  =  0.8-10“3  (3),  ^  =  10-3  (4), 
P  =  1.2-10  (5).  It  is  possible  to  remark  that  the  amplitudes  of  velocity  disturbances  are  small. 

Conclusions 

The  following  conclusions  can  be  drawn  from  the  present  investigations. 

1.  External  vortical  and  nonvortical  waves  entrained  by  the  stream  can  excite  high- 
intensity  disturbances  in  the  boundary  layer. 

2.  The  efficiency  of  disturbance  generation  increases  with  decreasing  frequency,  and  its 
maximum  is  observed  for  waves  with  a  characteristic  value  of  the  wavenumber  P*  in  the 

lateral  direction.  In  the  low-frequency  range,  the  value  of  p*  significantly  exceeds  the 
streamwise  wavenumber  a ,  which  is  characteristic  of  experimentally  observed  streamwise 
structures. 

3.  For  stationary  disturbances  and  M  =  0,  under  changes  in  Re  and  P  =  const ,  the  highest 
efficiency  of  disturbance  generation  in  the  boundary  layer  corresponds  to  the  value 

(PRe)i  «  0.7  ,  and  under  changes  in  p*  and  Re  =  const  corresponds  to  (PRe)2  «  0.55  .  The 
latter  value  is  in  good  agreement  with  the  period  of  streamwise  structures  obtained  in  the 
experiments  of  [11]. 

4.  The  efficiency  of  disturbance  generation  in  a  supersonic  boundary  layer  is  lower  than 
that  in  a  subsonic  boundary  layer. 

5.  Temperature  external  disturbances  excite  velocity  disturbances  inside  a  boundary  layer 
with  amplitudes  much  less  than  in  the  case  of  the  external  hydrodynamics  disturbances. 

The  authors  would  like  to  express  their  gratitude  to  G.V.  Petrov  for  the  calculation 
program  he  kindly  offered. 
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